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Fast radio bursts (FRBs) are extragalactic radio 
flashes of unknown physical origini] Their high lu- 
minosities and short durations require extreme energy 
densities, like those found in the vicinity of neutron 
stars and black holes. Studying the burst intensi- 
ties and polarimetric properties on a wide range of 
timescales, from milliseconds down to nanoseconds, is 
key to understanding the emission mechanism. How- 
ever, high-time-resolution studies of FRBs are limited 
by their unpredictable activity levels, available instru- 
mentation and temporal broadening in the interven- 
ing ionised medium. Here we show that the repeating 
FRB 20200120E4 can produce isolated shots of emis- 
sion as short as about 60 nanoseconds in duration, with 
brightness temperatures as high as 3 x 10*! K (excluding 
relativistic effects), comparable to ‘nano-shots’ from the 
Crab pulsar. Comparing both the range of timescales 
and luminosities, we find that FRB 20200120E observa- 
tionally bridges the gap between known Galactic young 
pulsars and magnetars, and the much more distant ex- 
tragalactic FRBs. This suggests a common magnetically 
powered emission mechanism spanning many orders of 
magnitude in timescale and luminosity. In this work, 
we probe a relatively unexplored region of the short- 
duration transient phase space; we highlight that there 
likely exists a population of ultra-fast radio transients 
at nanosecond to microsecond timescales, which current 
FRB searches are insensitive to. 

In the late 1960s, the discovery of radio pulsars was en- 
abled by observations using an unprecedented time resolu- 
tion of 50 — 100 mE! The subsequent discovery of the Crab 


pulsar®l eventually led to the insight that its giant pulses El 


show structure on timescales as short as 0.4 nanoseconds” 
These ‘nano-shots’ have exceptionally high brightness tem- 
peratures of ~ 107° KO The more recent discovery of the 
fast radio burst (FRB) phenomenon MMA has shown that na- 
ture produces millisecond duration radio flashes that are ap- 
parently 10°-!* times more luminous than those of average 
pulsars #3) 

As with pulsars, the recording of raw voltage data can 
allow us to study FRB signals on timescales of microsec- 
onds (us) down to nanoseconds (ns). Thus far, such ob- 
servations have demonstrated that FRB signals can show 
temporal structure on timescales of tens of microseconds 
down to just a few microseconds ES] Nimmo et al.“ also 
discuss how individual FRBs can display a range of tem- 
poral structure, from a few microseconds up to millisec- 
onds. The temporal behaviour of FRBs — along with their 
polarimetric properties 42] and dynamic spectra!8l — pro- 
vide key inputs for distinguishing between the dozens of 
proposed source models and emission mechanisms“. Fur- 
thermore, these burst properties can aid in observationally 
distinguishing the populations of repeating and apparently 
one-off FRBs, thereby a le the case that they have 
physically distinct origins 2222] 

Kirsten et al. (submitted) recently associated the repeat- 
ing FRB 20200120E4! with [PR95] 30244, a globular clus- 
ter that is part of the M81 galactic system. At a dis- 
tance of 3.6Mpc, FRB 20200120E is the closest-known 
extragalactic FRB, bridging the gap between the putative 
Galactic FRB source SGR 1935+2154 (a known magnetar 


at a distance of ~ 3 — 10kpc)2#25l and FRB 20180916B 
at a luminosity distance of 149 Mpc’, FRB 20200120E 
is also at high Galactic latitude (~ 41.2°), which reduces 
the effect of temporal scatter broadening arising from the 
Milky Way interstellar medium (ISM). This suggests that 
FRB 20200120E could be an excellent source to study at 
the highest-possible temporal resolutions. 

Here we present a spectro-polarimetric study of 5 
FRB 20200120E bursts detected with the Effelsberg 100-m 
telescope during a broader localisation campaign (Kirsten 
et al. submitted) using ad hoc interferometric observations 
with dishes from the European Very-long-baseline interfer- 
ometry (VLBI) Network (EVN). In this Article we present 
both extremely high time resolution, and lower time res- 
olution analyses of the bursts, with full polarimetry, and 
thereafter discuss the observational connections between 
FRB 20200120E and well-studied repeating FRBs, the Crab 
pulsar and the Galactic magnetar SGR 1935+2154. 


Observations & Data 


We observed FRB 20200120E“lusing an ad hoc EVN array 
at 1.4GHz, on 2021 February 20 UT 1700 — 2200, 2021 
March 7 UT 1545 — 2045, and 2021 April 28 UT 1100 - 
2200. For details on the interferometric array configura- 
tion and localisation results, see Kirsten et al. (submitted). 
With the 100-m Effelsberg telescope, we recorded dual cir- 
cular polarization raw voltages (R and L) using 32 MS/s real 
sampling per 16-MHz subband and 2-bit samples written in 
VDIF format=, with a total bandwidth of 256 MHz. 

Throughout this work we label the FRB 20200120E 
bursts as Bn, ordered according to their arrival time, 
and matching the nomenclature used in Kirsten et 
al. (submitted). 


Results 


As described in Methods, we create total-intensity filter- 
bank data containing each burst at the native sampling 
of the voltage data (31.25ns), using SFXC. The filter- 
bank data were created with 32-bit digitization, in order 
to avoid saturation. The data are coherently dedispersed 
within the 16 MHz subbands, and each subband is also time 
shifted to correct for dispersive delay (incoherent dedisper- 
sion), both using our measured dispersion measure (DM) of 
87.7527+0.0003 pc cm~> (Methods). We assume the same 
DM for all bursts. This DM is > 9c lower than the pre- 
viously reported measurement of 87.782 + 0.003 pe cm~? 
(where the quoted value is based on the average of 3 
bursts) The observed difference in DM could be due to 
unresolved time-frequency structure in the Bhardwaj et al. 
bursts, often seen in repeating FRBs!8l_ Future measure- 
ments are needed to determine if the DM is frequency or 
time dependent. 

In Figure[Ip, we present the profile of burst B3 at 31.25 ns 
resolution (black), and downsampled to 1 us resolution 
(green), in the frequency range 1318 — 1334 MHz. This 
range corresponds to the single subband containing the 
brightest spectral feature in the burst (visible in the burst 


dynamic spectrum shown in Figuref2}). By using a single 
subband, we avoid artefacts due to the inaccuracy of inco- 
herent dedispersion. There are clear few-bin-wide temporal 
structures in the 31.25 ns profile of burst B3 (Figure[Ih. The 
question remains whether the sub-microsecond structures 
we observe are isolated shots of emission or noise fluc- 
tuations consistent with the y7-distribution of amplitude- 
modulated noise (AMN; Methods). We calculate the prob- 
ability of drawing the bright ~60ns duration feature at 
Time = O in Figure[]] from the local brightness distribu- 
tion, where in this case “local” is defined as +1.5625 us 
(+50 bins) around the bright feature (Methods). We find 
that the probability of drawing this high signal-to-noise 
(S/N), two-bin-wide structure from the distribution is p= 
4x 10-8 x 100 bins/2 = 2 x 10~° (Extended Data FigurelI) 
Therefore, we find this feature to be inconsistent with the 
local AMN distribution, and conclude that this structure is 
a real isolated shot. In addition to this bright 60 ns shot, we 
find at least another 2 significant sub-microsecond shots of 
emission in burst B3. 

Bursts B2 and B4 also have sufficient S/N to study at the 
highest-possible temporal resolution (Methods; Extended 
Data Figure[]). Contrary to burst B3, which exhibits few- 
bin temporal structure, the highest S/N spikes in B2 and 
B4 are single-bin unresolved features. This poses a con- 
cern since the bright spectral feature dominating the sub- 
band used to create the burst profile, has a spectral extent 
less than the subband width (we attribute the spectral fea- 
tures to scintillation, see below). This results in an effective 
time resolution less than the native sampling of the data. 
Therefore, single-bin features are more likely to be consis- 
tent with a noise process. The probability of drawing the 
bright single-bin features from their local distributions are 
p=1x10-+x 100 bins = 0.01 and p = 4x 1074 x 100 bins = 
0.04 for B2 and B4, respectively. These relatively high 
probabilities, combined with the effective resolution argu- 
ment above, suggests that the high-resolution features in 
both B2 and B4 are consistent with the y? AMN distribu- 
tion. 

The total burst duration and spectral structure were 
quantified by performing a 2-dimensional autocorrelation 
of the lower-time-resolution burst dynamic spectra (8 us, 
125kHz; Methods). The burst temporal width and fre- 
quency extent are reported in Table[]] As is clear in the 
autocorrelation functions (ACFs; Extended Data Figure[2), 
there is an additional narrower frequency scale, which we 
measure to be consistent with the expected scintillation 
from the Milky Way ISM28l (Methods; Extended Data Fig- 
ure|3h). The scintillation bandwidth measurements are re- 
ported in Table[i] Additionally, we report the fluence, peak 
flux density and isotropic-equivalent spectral luminosity of 
the bursts, computed within the +20 width region (Table[I} 
Methods). 

There are three clear timescales measured in the aver- 
age temporal ACF (averaged over the four brightest sub- 
bands; Methods) of burst B3 at 31.25ns resolution (Ex- 
tended Data Figure/4f—h): a 40 us timescale, consistent with 
the total burst extent in time, a clear | us timescale, and 
even a shorter timescale (40 ns) consistent with temporal 
structure on the few-bin level. In contrast, bursts B2 and 
B4 both exhibit a timescale on the order of 10 us, con- 


sistent with their total burst duration (Extended Data Fig- 
ures [5] & (6). There is evidence for structure on the few- 
bin level in the B4 ACF, although the height of this narrow 
ACF feature relative to the wider ACF feature is smaller 
for B4 than B3 (height of cyan Lorentzian relative to green 
Lorentzian in Extended Data Figures 4h & [62), implying 
either that the S/N of these temporal fluctuations are lower, 
or that there are fewer temporal features on this timescale. 
No additional short-timescale components were measured 
in burst B2 (Extended Data Figure[5). This further supports 
the conclusions above that we have resolved microsecond 
and sub-microsecond structure in burst B3, with no evi- 
dence for similar structure in B2 and B4. Additionally, 
there are a range of timescales observed in the bursts, some- 
times observed within bursts: from tens of nanoseconds to 
tens of microseconds. In Figure[3}we compare the isotropic- 
equivalent luminosity of the three highest-significance shots 
of emission in the B3 31.25ns profile with other short- 
duration transients, including that of the bright ~ 5 us com- 
ponent of B3 seen in the 1 us resolution data (Extended 
Data Figure{7p), and the wider burst envelopes seen in B1, 
B2, B4 and B5 (Methods; Table[ī}. 


The power spectra of bursts B2, B3 and B4 (panel c; Ex- 
tended Data Figures[5}{6) are all consistent with red noise 
(Methods), with B2 and B4 exhibiting a steeper power law 
(a = 1.85 + 0.04 and 2.04 + 0.05, respectively) than B3 
(a = 1.46 + 0.05). This is consistent with the results we 
have presented above; the power spectrum of B3 shows 
more power at higher frequencies (shorter timescales), than 
bursts B2 and B4. 


As described in Methods, the polarimetric data were cal- 
ibrated using the known polarization properties of the pul- 
sar PSR B0355+5421 The full polarimetric burst profiles 
and polarisation position angles (PPA) are shown in Fig- 
ure[2] Note that we could not recover the polarimetric prop- 
erties of burst B5, likely due to the low S/N of the burst. 
Bursts B1 — B4 are highly linearly polarized (> 90%), and 
exhibit little-to-no circular polarization (Table2}. AS re- 
ported in Table] there is a tentative 3 — 4o detection of 
13 % and 6% circular polarization in B2 and B4, respec- 
tively. In the 8 us resolution profiles, there is evidence for 
small variations in the PPA across the bursts, with a APPA 
between bursts from the same epoch within ~ 30°. We de- 
termine the rotation measure (RM) of the bursts, and con- 
clude that they are in agreement with previous measure- 
ments (—29.8 rad m24 


As presented above, burst B3 exhibits variations on both 
microsecond as well as sub-microsecond timescales. The 
Stokes parameters are only physically meaningful with suf- 
ficient averaging? Therefore, in Extended Data Figure|8| 
we show the frequency-averaged polarization profile and 
PPA of burst B3 at 125ns resolution. In these data we av- 
erage over a total of 44 channels (over the frequency range 
1254 — 1430 MHz), and therefore have 44 degrees of free- 
dom. The ~100 ns — 1 us structures are highly linearly po- 
larized, consistent with the polarization properties at lower 
time resolution. Additionally, the PPA varies between the 
sub-us temporal features by up to a few 10s of degrees. 


Discussion 


The timescales and luminosities measured in the 
FRB 20200120E bursts presented in this work popu- 
late a previously vacant, relatively unexplored region of 
the short-duration transient phase space. They bridge the 
gap between extragalactic FRBs and Galactic neutron stars 
(Figure[3). Here we elaborate on how FRB 20200120E 
compares observationally with other short-duration radio 
transients, and discuss the implications of our findings. 


Repeating fast radio bursts 


The polarimetric properties of FRB 20200120E are consis- 
tent with those of most well-studied repeaters 46117] As is 
often observed, we find a very high (~ 100%) linear and 
low (<$ 10%) circular polarization fraction as well as a 
polarization position angle (PPA) that is roughly constant 
during bursts. As with FRB 20180916B, we find that 
FRB 20200120E shows small PPA variations (APPA) on 
timescales < 10 us. Between bursts, we see variations in 
the PPA of a few tens of degrees, in contrast to the < 10° 
APPA from bursts detected at the same observing epoch, 
seen in other repeaters 4116) We note that at least one re- 
peater (FRB 20180301A) has shown lower linear polariza- 
tion fractions (40 — 80%) and significant PPA swings in 
some bursts] 

The spectrum of FRB 20200120E shows at least two 
scales of variation. We find narrow-band brightness varia- 
tions on the scale of ~ 5MHz, which we ascribe to scin- 
tillation in the Milky Way ISM. The other ~ 100 MHz 
spectral variation may be intrinsic to the source emission 
mechanism or imparted by local propagation effects. This 
spectral envelope is similar to what is seen from other re- 
peaters i8222 Additionally, the downward-drifting burst 
sub-structure, often referred to as the ‘sad-trombone effect’, 
that is often seen in repeaters, has previously been ob- 
served for FRB 20200120E. 

Though FRB 20200120E shares many characteristics of 
repeating FRBs, its ~ 100 us burst envelopes are atyp- 
ically narrow. The range of timescales observed, how- 
ever, roughly a factor of 1000, is comparable to what has 
been found in a similar analysis of FRB 20180916B% 
In the case of FRB 20180916B, Nimmo et al.[6l were 
limited by a larger scatter-broadening of 2.7 us (note that 
FRB 20180916B is at a Galactic latitude of only 3.7°; 
see also Extended Data Figure Bhe and could not rule 
out the possibility that the wider sub-bursts are composed 
of closely spaced microsecond structures. Burst B3 from 
FRB 20200120E has clear isolated shots of duration ~ 
60 ns, and also evidence that the sub-microsecond shots are 
clustered on microsecond timescales (Extended Data Fig- 
ure[4f-h). For bursts B2 and B4, there is no clear evidence 
for isolated shots in the 31.25 ns burst profile, but it is pos- 
sible that the S/N is too low to detect these individual shots, 
consistent with the lower <0.33 measured correlation coef- 
ficient (Methods; Extended Data Figure{7}i-). 

The strikingly similar observational properties of 
FRB 20200120E with other repeaters suggests that they 
have similar physical origins. FRB 20200120E is 
the closest known extragalactic FRB discovered to date 


(Kirsten et al. submitted). The remarkable proximity of 
FRB 20200120E has revealed radio bursts with an isotropic- 
equivalent spectral luminosity 2-3 orders of magnitude 
weaker than bursts from other repeating FRB (Figure[3). 
Such low luminosity bursts would not be detectable at the 
distance of any other precisely localized repeating FRBs. 
Continued monitoring of FRB 20200120E will be impor- 
tant to compare the energy distribution and activity rate with 
other repeaters. 


Crab pulsar 


We have discovered resolved structure in bursts from 
FRB 20200120E with durations of ~ 5 us, down to ~ 60 ns, 
2 orders of magnitude shorter timescales than have been 
probed for FRBs, to date. Crab pulsar giant pulses (GPs) 
show temporal structure in the range us — ms, which some- 
times resolve down to sub-nanosecond timescales, often re- 
ferred to as ‘nano-shots’2 Crab nano-shots are known to 
be extremely energetic, with brightness temperatures re- 
ported up to 10*! K (ignoring relativistic effects) The 
sub-microsecond shots in B3 exhibit a comparable, ex- 
tremely high brightness temperature (again, ignoring rela- 
tivistic effects; Figure B). This is at least four orders of 
magnitude higher than typically seen from the FRB popula- 
tion (10°2-1037 K)4 In the line of sight to FRB 20200120E, 
the expected scatter broadening from the Milky Way ISM 
is ~ 50ns at 1.4GHz28) in rough agreement with our mea- 
sured scintillation bandwidth (1/(27Ay¢cint) ~ 27 ns). This 
implies that the ~ 60ns structure observed in this work is 
likely the shortest resolvable temporal scale at our observ- 
ing frequency. 

Bursts B2, B3 and B4 are consistent with the Scintil- 
lating Amplitude Modulated Polarized Shot Noise model 
(Methods), which has been used to describe many aspects 
of pulsar emission, including GPs from the Crab pulsar3324] 
However, only B3 is observed to resolve down to microsec- 
ond and sub-microsecond temporal scales. This is simi- 
lar to observations of the Crab pulsar, where not all GPs 
resolve down to individual nano-shots/235) Furthermore, 
most Crab GPs consist of at least one broadband ‘micro- 
burst’, with a characteristic timescale of a few microsec- 
onds0361 similar to the structure that can be seen in burst 
B3 at 1 us time resolution and consistent with the measured 
1.11 us timescale in the temporal ACF (Extended Data Fig- 
ures|4b &[7p). Sometimes, the Crab micro-bursts are seen 
to resolve down to individual narrowband nano-shots (see 
e.g. Figure 4 of Hankins & Eilek#%), consistent with the 
features seen at 31.25 ns in burst B3. 

Band-limited giant pulses from the Crab pulsa: [BBS] as 
well as from the ‘Crab’s twin’ pulsar PSR J0540-691922) 
have been observed. These are reminiscent of the narrow- 
banded emission of observed from repeating FRBs. This 
adds further weight to the phenomenological connection of 
giant pulse emission with FRBs. 


Galactic magnetar SGR 1935+2154 


The FRB 20200120E bursts also exhibit wider charac- 
teristic timescales of 10 — 100 us. These wider compo- 
nents have isotropic-equivalent energies on the order of 


10°2-33 erg, which is 2-3 orders of magnitude weaker than 


the bright FRB-like radio burst from the Galactic magne- 
tar SGR 1935+21542324 Thus, SGR 1935+2154 has pro- 
duced more energetic radio bursts than some extragalac- 
tic FRBs, eliminating the gap in luminosity (Figure[3) and 
thereby strengthening the connection between magnetars 
and FRBs. The observed energies are much lower than the 
proposed low-energy cutoff at ~10°* erg for the low-twist 
magnetar model) implying that the FRB luminosity func- 
tion does not abruptly end at ~10** erg. The brightness tem- 
peratures of these wider components are ~ 10°? K, com- 
parable to SGR 1935+2154 231241 and also consistent with 
the lower end of the observed FRB brightness temperatures 


(Figure B}. 


Implications 


Previously, it has been shown that millisecond-duration ra- 
dio bursts from the Galactic magnetar SGR 1935+2154 
span 7-8 orders of magnitude in apparent luminosity ®, 
bridging from the brightest pulsars, to 1—2 orders of magni- 
tude lower than the weakest known extragalactic FRBs. By 
probing FRB 20200120E at timescales of tens of nanosec- 
onds, we highlight that this coherent radio source observa- 
tionally ties repeating FRBs with GP emission from young 
pulsars, and the bright FRB-like emission from magnetars. 
The timescales and luminosities measured in this work fill 
the gap in the luminosity-duration phase space (Figure[3), 
further emphasising that observationally the division be- 
tween source populations (FRBs, pulsars, magnetars) is not 
clear. 

Constraints on the shortest timescale variations in FRB 
lightcurves are key for understanding the physical mecha- 
nism producing the bursts, and can ultimately reveal clues 
about the progenitor. The observational connection to the 
Crab pulsar, SGR 1935+2154, and repeating FRBs, sup- 
ports a common magnetically-powered emission mecha- 
nism spanning many orders of magnitude in timescale and 
luminosity. The ~ 60 ns to 5 us timescales observed in B3 
correspond to a light-travel size of 20-1500 m, ignoring rel- 
ativistic effects. The sub-microsecond timescales observed 
are too short to be naturally explainable via emission from 
a synchrotron maser in a relativistic shock22) since it would 
require a small region to be emitting at any given time. 
Nimmo et al. previously argued that the short timescales 
observed in FRB 20180916B are more naturally explain- 
able in the context of a neutron star magnetospheric origin. 
The results presented here further support a magnetospheric 
origin of the FRB emission. 

Magnetic reconnection is one possible way to harness the 
magnetic energy to power a coherent radiation mechanism 
that produces the radio emission. This has been proposed to 
explain Crab nano-shots=2) the FRB-like radio burst from 
SGR 1935+21544) and FRBSE] A range of temporal 
scales and magnetic energy releases can be expected from 
this mechanism, which is consistent with the observed dy- 
namic range of timescales, and large range of FRB lumi- 
nosities; from the weaker possible FRBs emitted by the rel- 
atively old magnetar SGR 1935+2154, and the bursts from 
FRB 20200120E, to the more energetic FRBs potentially 
coming from extragalactic young, active magnetars. 


The observed timescales and luminosities from 
FRB 20200120E can be explained by magnetic re- 
connection events in the close vicinity of a relatively 
young, isolated, highly magnetized neutron star. The 
association of FRB 20200120E with an old globular cluster 
strongly implies that, if the source is a magnetar, it was not 
formed through a core-collapse supernova (Kirsten et al. 
submitted). The globular cluster origin of FRB 20200120E 
also allows for the exploration of alternatives to a magnetar 
progenitor: for example, an accreting system, in which 
case the observational similarities with the Crab pulsar and 
SGR 1935+2154 are more coincidental. 

Future observations of FRB 20200120E at observing 
frequency > 5GHz (where the scatter broadening will 
be lower), with a bandwidth > 200MHz, are needed 
to explore shorter timescales. Continued monitoring of 
FRB 20200120E will provide statistics on the distribution 
of emission timescales and whether the burst activity is pe- 
riodic, like FRB 20180916B22 Such data will also help 
determine whether, e.g., there are quasi-periodic fluctu- 
ations in the burst lightcurves, hinted at in the case of 
FRB 20180916826] or a secular variation in the widths of 
burst envelopes. 

Lastly, we highlight that the short timescales measured 
from FRB 20200120E in this work motivate searches for a 
population of ultra-fast radio bursts, despite the consider- 
able technical challenges that such a search entails. 
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Figure 1: Burst B3 from FRB 20200120E exhibits sub-microsecond temporal structure. Panel a shows the temporal profile 
of burst B3 at 31.25 ns time resolution (black) and downsampled to 1 us resolution (green). Panel b shows a histogram of 
the S/N values in the 31.25 ns data (black), also showing the x-distributed (2 degrees of freedom) noise distribution (red) 
arising in the limiting case of ôtôv ~ 1, where ôt and ôv are the time and frequency resolutions of the data, respectively. 
The residuals of the off-burst histogram with the best-fit y?-distributed with 2 degrees of freedom is shown in panel c. 
This plot contains temporal profiles computed for the bright scintle only, i.e. the frequency range 1318 — 1334 MHz. The 
data were generated with SFXC and are coherently dedispersed to 87.7527 pe cm™°?. Panel d shows a zoom in of panel a, 
highlighting the temporal regions plotted in panels e (yellow) and f (blue) panels. 
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Figure 2: The polarimetric profiles, dynamic spectra, time-averaged spectra and polarization position angle (PPA) of 
the bursts detected from FRB 20200120E. For B5 only Stokes I is shown (see text). The data are plotted with 8 us 
and 1 MHz time and frequency resolution, respectively (with the exception of B1 and B5 which are plotted with 4 MHz 
frequency resolution). The data were generated with SFXC and are coherently dedispersed within each 16 MHz subband to 
87.75 pe cm~? (and also incoherently shifted between subbands). Panels a-d are the PPA across the burst profile, where the 
color gradient represents the linear polarization S/N (black is high S/N and white is low S/N), and the colored horizontal 
line represents the weighted best-fit line to the PPA. Only the PPAs above a linear S/N threshold of 3 are plotted. Panels 
e-h show the total intensity (Stokes I; black), unbiased linear polarization (Everett & Weisberg8h red) and circular 
polarization (blue) burst profiles (panel i shows the Stokes I profile of burst B5). In the top-left of the panels are the burst 
name used throughout this work, and the time resolution used for plotting. The colored bar at the bottom of the panels 
represent the +20 burst width used to measure the polarization fractions and burst fluence. Panels j-n are the dynamic 
spectra and panels o-s are the time-averaged frequency spectra. The red marks on the dynamic spectra outline the edges 
of the subbands. Data that have been removed due to radio frequency interference have not been plotted. 
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Figure 3: Nanosecond to second transient phase space. The FRB 20200120E burst temporal structures and their measured 
isotropic-equivalent spectral luminosity are shown by the black stars. We plot the three highest significance isolated 
shots from burst B3 (left-most black stars), the bright 5 us structure evident in the 1 us time resolution data of burst 
B3 (Extended Data Figure), as well as the wider burst structure seen in B1, B2, B4 and B5 (right-most black stars). 
The other FRBs plotted are the published localized repeating FRBs, with known distances: FRB 20121 102A 1408625158) 
FRB 201 809168 SZU2671960 and FRB 2019071 ARIE In purple we plot the radio bursts from the Galactic magnetar 
SGR 1935+2154232441l63l64) Giant pulses from the Crab pulsar (GRPs) are shown in orange and the ‘nano-shots’ 
are represented by the yellow crosses EBS] The pulsar and Rotating RAdio Transient (RRAT) population are shown in 
pink®*! Note that for pulsars, the points on the figure represent average pulses, where individual pulses can be brighter or 
weaker by approximately 1 or 2 orders of magnitude. The grey lines represent lines of constant brightness temperature. 


Table 1: Burst properties. 


Burst Time of Arrival* Fluence?” Peak S/N° Peak Flux Density®* Spectral Luminosity>** Width T Frequency Extent’ Scintillation bandwidth" 
[MJD] [Jy m0s] [Jy] [10% erg s7! Hz7!] [us] [MHz] [MHz] 

Bl 59265.88304437179 0.13 + 0.03 6.6 1.59 + 0.32 0.56 + 0.13 156+ 1 140 +1 1.9 +0.7 

B2 59265.88600912486 0.63 +0.13 36.1 8.71 + 1.74 2.3+0.6 62+1,93+0.58 103+ 1,89+18 3.4+ 1.0 

B3 59280.69618745651 0.53+0.11 64.8 15.6 + 3.12 8.0+1.9 46.7 + 0.1 9441 6.9 + 1.1 

B4 59280.80173397988  0.71+0.14 29.3 7.07 + 1.41 4.0+1.0 117+1 134 +1 6.8 + 1.3 

BS 59332.50446581106 0.09 + 0.02 6.9 1.66 + 0.33 1.0+0.2 56.6+0.1 8641 29415 


a Corrected to the Solar System Barycenter to infinite frequency assuming a dispersion measure of 87.75 pe cm7, 


reference frequency 1502 MHz and dispersion constant of 1/(2.41x10-+) MHz? pe~! cm? s. 
The times quoted are dynamical times (TDB). 
> We estimate a conservative 20% error on these measurements, arising due to the uncertainty in the system equivalent flux density (SEFD) of Effelsberg. 
€ Determined for a time resolution of 8 us. 
4 Using the distance to the globular cluster [PR95] 30244 of 3.63 +0.34 Mpd®&] 
e Isotropic-equivalent. 
f Defined as 1/ V2 multiplied by the FWHM 
of the autocorrelation function (ACF). 
£ Width per burst component. 
The uncertainties are a combination of the 1-ø fit uncertainty, and 1/ WN, where N ~ 5 is the number of scintles. 


Table 2: Burst polarization properties and polarization position angle fit results. 


Burst RM? [radm?] PPA offset? [deg] RM (global)** [rad m7] PPA offset (global)° [deg] Lunbiasea/I9[%] V/IT[%] y7° Degrees of freedom? 
Bl —21.9 + 13.1 28.5 —54.2 + 4.8 68.3 9449 4+9 18.9 T 
B2 -57.2 + 5.1 0 —54.2 + 4.8 0 103 + 3 -13+3 596 27 
B3 —37.1 + 4.2 —8.6 —36.9 + 3.8 -9.3 102 +2 1+1 15.5 11 
B4 —36.4 + 8.0 13.1 —36.9 + 3.8 14.3 9442 6+2 25.7 22 


* The large fractional errors arise due to covariances between fit parameters in the QU-fitting, see text for more details. 


> Offset from a weighted mean of the bursts from a single epoch. 


€ The global values are assuming all bursts detected on the same day have the same RM. 
4 The quoted uncertainties are 1a uncertainties, and do not include uncertainties from the calibration. 
Note: we reproduce the published pulsar circular (linear) fractions within 3% (10%) and 1% (6%) 


for the Feb 20 and March 7 observations, respectively. 
€ Fit results of weighted least squares fitting of a straight line to the PPA. 


Methods 


In this section, we describe the analysis performed in more 
detail. Since this was an interferometric campaign, we reg- 
ularly observed a phase calibrator source (5.5 minutes on 
target, 2.0 minutes on the calibrator). The total time on 
FRB 20200120E was 2.93 hr per observation. 


Burst search 


We performed two independent analyses on the single-dish 
Effelsberg data to search for bursts from FRB 20200120E. 
The voltage data were converted to total-intensity filter- 
banks using digifil with time and frequency resolu- 
tions of 64 us and 125 kHz, respectively. These data were 
searched using a Heimdall-based pipeline and burst can- 
didates were classified using the machine learning clas- 
sifier, FETCH] A more detailed description of this 
pipeline was presented in Kirsten et ale In parallel, we 
also recorded total intensity filterbanks using the PSRIX 
pulsar backend®8l with time and frequency resolution of 
102.4 us and 0.49 MHz, respectively, and a total bandwidth 
of 250MHz. The PSRIX pulsar backend data were anal- 
ysed using a PRESTO-based pipeline) and single pulse 
candidates were classified using an automated classifier 
based on Michilli et al We recorded a test pulsar scan of 
PSR B0355+54 to inspect the data quality. 


Five bursts from FRB 20200120E were discovered in 
the search of the raw voltages using the Heimdal1-based 
search (B1 and B2 on 2021 February 20; B3 and B4 on 
2021 March 7; and B5 on 2021 April 28). Only four of 
the five bursts were found in the PSRIX data because burst 
B1 occurred during the ~ 40s latency between the start of 
the VLBI recording and the pulsar backend recording. Ini- 
tially, B2 was flagged as radio frequency interference in the 
PSRIX data by the automated classifier due to strong scin- 
tillation structure. 


Data products 


For further analysis of the burst properties — using full po- 
larimetric information at a range of time and frequency res- 
olutions — we generated filterbank data from the raw volt- 
ages using both SFXC“ and digifil™| 


SFXC 


Using SFXC (phased array branct{'), we produced rela- 
tively low time resolution filterbank data (8 us), in order 
to have sufficient spectral resolution (125 kHz) to study fre- 
quency structure due to scintillation, and to determine the 
rotation measure (RM). The data are coherently dedispersed 
within the 16-MHz subbands, and incoherently dedispersed 
between the subbands using a DM of 87.75 pecm™, which 
is within ~ 0.003 pc cm~? of the value we measured using 
higher time resolution data (see below). This 0.003 pe cm~? 
offset corresponds to a dispersive delay, across the burst ex- 
tent in frequency, that is less than the time resolution. Using 
dspsr2l we created archive files containing each burst at 
the resolution of the filterbank data, and remove frequency 
channels from the data that are contaminated by radio fre- 
quency interference using the psrchive tool pazil 

The real-sampled VDIF data with 16-MHz subbands al- 
low for up to 31.25 ns time resolution. Due to the Fast 
Fourier Transform-based (FFT-based) correlation and de- 
lay corrections implemented in SFXC (phased array), the 
practical time resolution limit is 125 ns. In this case, the 
data have 4 channels per subband. Since SFXC does not 
window the data before the FFT, potential spectral leakage 
artefacts are a concern, especially in the case of a low num- 
ber of channels per subband. This data is used to study the 
polarimetry at high time resolution. 

Additionally, we output coherently dedispersed total in- 
tensity (Stokes I) at the original data sampling rate (i.e. no 
channelization within subbands) using the SFXC bolome- 
ter brancl}?| Here we have no concerns of spectral leak- 
age or any other artefacts that may be introduced at the 
channelization stage of the processing. Note that the VLBI 
backed, DBBC24 uses digital down conversion to cre- 
ate independent 16 MHz subbands, meaning that adjacent 
31.25 ns samples are not correlated. 


Digifil 


We also used digifil to produce data at the highest- 
possible time resolution, in order to verify the results pre- 
sented in this work. digifil utilises a convolving filter- 
bank, which reduces spectral leakage effects 2 However, 
digifil converts to complex-sampled data internally, re- 
sulting in a factor of 2 decrease in the time resolution for a 
given number of channels (for the Nyquist limit). Thus, the 


l https://github.com/aardk/sfxc/tree/phased-array 
?https: //github.com/aardk/sfxc/tree/bolometer 


best time resolution we can achieve with our 16-MHz sub- 
bands is 62.5ns. digifil produces one filterbank file per 
subband, and the subbands are combined with incoherent 
dedispersion. 

In both the digifil and SFXC data, the incoherent 
dedispersion between subbands is performed with integer 
sample shifts. 


Refined dispersion measure determination 


At the native time resolution of the voltage data (31.25 ns), 
the DM must be correct to < 0.0005 pe cm™ to avoid DM 
smearing across the 256 MHz bandwidth larger than the 
time bin width. 

To limit potential spectral leakage artefacts in the SFXC 
data, we produced 500 ns time resolution filterbank files, 
with 16 channels per subband (1 MHz channels), at a range 
of DM values around the expected DM (determined by eye 
in the dynamic spectrum of burst B3). The data are coher- 
ently dedispersed within the 16 MHz subbands, and each 
subband is then time shifted to correct for dispersive de- 
lay, to limit DM smearing within each channel, and across 
the band. These data were then downsampled in time to 
1 us (i.e. a downsampling factor of 2). The reason for pro- 
ducing the data products in this way, as opposed to cre- 
ating 1 us resolution data with higher frequency resolu- 
tion, is to minimise smearing due to incoherent dedisper- 
sion between subbands (integer bin shift). We can accu- 
rately determine the DM of burst B3 by maximizing the 
S/N of the bright ~ 5 us component (Extended Data Fig- 
ure[9), despite any downward-drifting, ‘sad-trombone’ ef- 
fect, that may be present in the burst dynamic spectrum. 
This is because by maximizing the S/N of the narrow tem- 
poral structure, we are essentially maximizing the structure 
in the burst (as opposed to maximizing the S/N of the en- 
tire burst, which has been shown to incorrectly capture the 
DM, e.g. Gajjar et al.2h, In Extended Data Figure p] 
we plot the peak S/N of the burst profile as a function of 
DM, and fit a Gaussian function to determine the best DM. 
We search for a DM in the range 87.74 — 87.765 pecm™> 
in steps of 0.001 pccm™ and find that the best DM is 
DM= 87.7527 + 0.0003 pecm™, where the uncertainty is 
determined by o//A, where o and A are the standard devi- 
ation and amplitude of the Gaussian fit, respectively. We 
note that the range of DMs searched do not include the 
previously measured DM since the dynamic spectrum of 
burst B3 is visibly over-corrected when dedispersed to this 
value. Both measurements of the DM use the same disper- 
sion constant of 1/(2.41x10-+) MHz? pe! cm? s, thus they 
can be directly compared. 

In panels c and d of Extended Data Figure p] we plot 
the burst profile and dynamic spectrum coherently (within 
channels) and incoherently (between subbands) dedispersed 
to this best-fit DM, respectively. This analysis was re- 
peated on digifil filterbank data products generated 
with 62.5ns time and 16 MHz frequency resolution, co- 
herently and incoherently dedispersed, and then down- 
sampled to | us time resolution. We find a consistent 
DM value from digifil and SFXC data (noting also 
that both softwares use the same dispersion constant of 
1/(2.41x107*) MHz? pe™! cm? s). 


We assume the same DM for all bursts in our sample. 


Scintillating Amplitude Modulated Polarized 
Shot Noise 


In the uncertainty principle limit ótôv ~ 1, where ôt and ôv 
are the time and frequency resolution (respectively) of indi- 
vidual samples, the off-burst noise is y°-distributed with 2 
degrees of freedom (Figure [Ip). For a modestly broadband 
noise-like signal, the statistics are y7-distributed where the 
degrees of freedom depend on the polarization fraction. We 
perform a least-squares fit of a y?-distribution to the on- 
burst distribution and find the best fit to have 2 degrees 
of freedom, consistent with 100% polarized Scintillating 
Amplitude Modulated Polarized Shot Noise (SAMPSN)2>! 
The superposition of the off-burst 77-distribution and the 
on-burst y7-distribution reasonably describes the total burst 
S/N distribution, where the residuals can likely be attributed 
to the fact that the shape of the burst envelope is not well- 
modelled. 

Due to strong spectral dips near the subband edges, suffi- 
cient frequency resolution is required in order to robustly 
correct the bandpass and compare individual bin spectra 
with each other. We therefore created a 1 us, 500 kHz fil- 
terbank using SFXC. Shown in Extended Data Figure [7]is 
the correlation coefficient between individual 1 us time bins 
above a S/N threshold of 9 in burst B3, as a function of their 
time separation. We find that the correlation coefficients 
have a geometric mean S/N weighted average of 0.31. The 
scintillation frequency structure is expected to be perfectly 
correlated within the duration of the burst (we observe that 
the B1 and B2 spectra, separated by 4.3 minutes, are corre- 
lated; Extended Data Fi gure[1 0p, while the Amplitude Mod- 
ulated Polarized Shot Noise frequency structure will change 
depending on the degree of polarization®4] SAMPSN pre- 
dicts a correlation coefficient of 


7 1 
~ 24 p2 


frac 


P (1) 
where Prac is the total polarization fraction. Our measured 
p = 0.31 implies that the signal is 100 % polarized, consis- 
tent with the high linear polarization fraction we measure 
in the frequency-averaged burst profiles (Figure B). and the 
on-burst S/N distribution. 

For B2 and B4, we find a weighted average correla- 
tion coefficient of the individual 1 us time bin spectra 
(above a S/N of 9) of 0.19+0.001 and 0.16 0.001, re- 
spectively, significantly lower than the 0.33 expectation for 
SAMPSNEŽ, and the 0.31 measured from B3, but also sig- 
nificantly greater than 0. This is consistent with highly 
polarized SAMPSN (similar to B3), where potentially the 
lower correlation coefficient can be attributed to low S/N 
and/or sparseness of shot pulses. This is consistent with the 
lack of evidence for resolved shots of emission in the high 
time resolution profiles of bursts B2 and B4 (Extended Data 


Figure[i}. 
Resolved sub-microsecond emission? 


To test whether the sub-microsecond temporal structures we 
observe are isolated shots of emission or, alternatively, con- 


sistent with AMN, we compare the high S/N features with 
the local brightness distribution. For burst B3, which ex- 
hibits envelope fluctuations on the us level (Extended Data 
Figure{7p), we define the local distribution as +1.5625 us 
(+50 bins) around the brightest feature (which we wish to 
measure the significance of). We take this range of bins 
to ensure that we are including enough samples to measure 
the local distribution, and to ensure we are not taking too 
many samples such that we lose the information about the 
burst envelope locally. Using an Anderson-Darling test 24) 
we confirm that this local distribution is exponentially dis- 
tributed. Note that since we cannot distinguish (by eye) be- 
tween temporal spikes that are due to AMN versus individ- 
ual emission spikes, we include all time bins (within the 
range mentioned), excluding the central feature, in the dis- 
tribution; therefore, the probabilities measured are lower- 
limits. We fit a y7-distribution with 2 degrees of free- 
dom, using a least-squares fit, to the S/N values within the 
range defining the local distribution. This distribution is 
expected for 100 % polarized SAMPSN|_ For the high 
S/N, 2-bin wide feature in burst B3 (Time= 0, Figure[Th, 
the probability of drawing this from the local distribution 
isp= 4x 1078 x 100 bins/2 = 2 x 10~° (Extended Data 
Figure[]), i.e. inconsistent with AMN, supporting that this 
component is a resolved, isolated shot of emission. 

We repeated the same analysis on bursts B2 and B4, since 
they have sufficient S/N to study at the highest time resolu- 
tion. For B2 we created the temporal profile at 31.25 ns 
using the subband 1302-1318 MHz, and for B4 we use sub- 
band 1398-1414 MHz: in both cases this is the subband 
containing the brightest spectral feature (Extended Data 
Figure[I). In the case of B3, at the highest time resolu- 
tion there are temporal fluctuations that are multiple bins 
wide. In contrast, the only bright structure in both B2 and 
B4 at the highest time resolution are unresolved single-bin 
spikes. Given that the frequency extent of the bright spec- 
tral feature is less than the subband width of 16 MHz (we at- 
tribute the spectral features to scintillation; see below), this 
results in an effective time resolution lower than the sam- 
pling resolution. Single bin spikes in this case are there- 
fore more likely to be consistent with the noise process. 
We tested the significance of the brightest unresolved struc- 
ture using the same method as described for B3. Again, 
we confirm that the local distribution around the structure 
of interest is exponentially distributed using an Anderson- 
Darling testé, The probability of the highest unresolved 
spike in B2 and B4 is p = 1 x 1074 x 100 bins = 0.01 and 
p= 4x 1074 x 100 bins = 0.04, respectively. Given the ef- 
fective resolution argument above, combined with the high 
probabilities of occurring by chance, there is not strong ev- 
idence supporting that these structures are isolated, unre- 
solved shots of emission. Instead, we conclude that these 
features are consistent with the y? AMN distribution. 


Temporal ACF and power spectrum 


In Extended Data Figures [4}{6] we present each coherently 
dedispersed subband of burst B3, B2 and B4, respectively, 
(covering the frequency extent of the burst) at 31.25 ns res- 
olution. This is essentially the burst dynamic spectrum at 
poor frequency resolution. We select the four brightest sub- 


bands that contain the most burst structure, and computed 
the ACF of the time profile. These ACFs were averaged 
together and shown in panel e of Extended Data Figures} 
[6] The reason for computing the ACF in this manner, as 
opposed to creating a frequency-averaged profile and com- 
puting the ACF, is to limit the smearing due to inaccuracies 
in the incoherent dedispersion (the subbands are shifted by 
an integer number of bins). For all three bursts, we measure 
a characteristic timescale on the order of 10 us, determined 
by fitting a Lorentzian function to the ACF (out to a time lag 
that is chosen by eye, to help the fitting distinguish between 
multiple timescales). This timescale is consistent with the 
full burst extent in time. In burst B3, we measure an addi- 
tional two timescales: a clear 1.11 us timescale, and even a 
shorter timescale (40 ns) consistent with temporal structure 
on the few-bin level (Extended Data Figure[4f-h). For B4, 
we additionally measure a short-temporal scale suggesting 
that there is temporal structure of a few-bins in the B4 
31.25 ns profile (Extended Data Figure|6p). Although, it is 
worth noting that the height of this narrow ACF feature rela- 
tive to the wider ACF feature (height of the cyan Lorentzian 
relative to the green Lorentzian in Extended Data Figures|4h 
and (6) is smaller for B4 than B3, implying either that the 
S/N of these temporal fluctuations are lower, or that there 
are fewer temporal features on this timescale. In B2, we 
see no evidence for power on shorter timescales in the ACF 
(Extended Data Figure[5p-f). 

The burst power spectra are presented in panel c of Ex- 
tended Data Fi gures|4}16] We perform a Bayesian maximum 
likelihood fit of a power law (red noise) plus constant (white 
noise) of the form 


f(y) =Av%+C, (2) 
where A is the amplitude, œ is the slope of the power law, 
and C is a white noise component, to the power spectrum, 
using the Stingray modelling interface, Additionally, we 
fit a red noise plus white noise plus Lorentzian model to 
the data, to search for the presence of a quasi-periodic os- 
cillation (QPO). This search was motivated by the hint of 
structure in the ACF of bursts from FRB 20180916816, and 
additionally, the low amplitude, wide frequency bump seen 
by eye in the downsampled power spectrum of burst B3 
(Extended Data Figurefþ). First, we compute the Bayesian 
Information Criterion for each model fit, defined as 


BIC = —2ln(L) + kln(n), (3) 
where L is the maximum likelihood of the fit, k is the num- 
ber of parameters in the model, and n is the number of data 
points. A lower BIC implies the data is better represented 
by that model, although it does not mean that the preferred 
model is a good fit to the data. Since the number of model 
parameters k is included in the BIC, models with more pa- 
rameters are penalized to avoid overfitting the data. Here 
we computed the A BIC = BICp;—-BICpr,),,, where BIC p;, 
is the BIC of the red noise model, and BICpz+or is the BIC 
of the red noise plus QPO model. For B3, we measure a 
ABIC= —12.1, which is significantly in favor of the red 
noise model. We then measure the goodness-of-fit p-value 
of the red noise model, which is the fraction of 100 simula- 
tions (using the MCMC package emcee) with a maximum 


10 


likelihood lower than the likelihood of our fit. This p-value 
is 0.53, implying the fit to the data is good. Finally, we test 
the significance of the highest outlier in the residuals of our 
best-fit model. The residuals are defined as 


2P(y) 


R) = Moy’ 


(4) 


for the power spectrum P(v) and best-fit model M(v). By 
simulating 100 power spectra from the posterior distribu- 
tion (using emcee) and computing the highest outlier for 
each simulation, we can compute a p-value (the fraction of 
simulated outliers that are higher than our measured out- 
lier). We find no significant outliers in the residuals, imply- 
ing we have no evidence of a QPO in the data. The power 
spectra of B2 and B4 are also found to be consistent with 
red noise. A summary of the A BIC, red noise goodness-of- 
fit p-value with power law index, and p-value of outliers are 
presented in Extended Data Table{I] 


Fluence and luminosity 


The burst profiles, in S/N units, are converted to physical 
units (flux density, Jy) using the radiometer equation” us- 
ing typical values for Effelsberg’s 1.4 GHz receiver temper- 
ature (20 K) and gain (1.54JyK~'). We expect these sys- 
tem values to be accurate to within 20%, which dominates 
the errors on the peak flux density and fluence (note that 
the peak flux density depends strongly on the time resolu- 
tion used, and the values reported in Table[i] are calculated 
using a time resolution of 8 us). We additionally consider 
a sky background temperature of 0.8 K, by extrapolating 
from the 408 MHz map, using a spectral index of -2.781 
and adding a 3 K contribution from the cosmic microwave 
background®4] For the 8 us burst profiles, we report the 
peak S/N, peak flux density and fluence (measured in the 
+2o width region) in Table[I] We also report the isotropic- 
equivalent spectral luminosity of the bursts, taking the dis- 
tance to FRB 20200120E as 3.63 + 0.34 Mpc (Kirsten et al. 
submitted). 


Burst temporal extent and spectral structure 


We performed a 2-dimensional Gaussian fit to the burst dy- 
namic spectra (Figure|2), to determine the burst extents in 
time and frequency. The measured widths were clearly un- 
derestimated compared with what can be seen by eye in 
the dynamic spectra, likely since the bursts are not well- 
modelled by a 2-dimensional Gaussian function. In the case 
of B1, B3, B4, and B5 we use the Gaussian mean in time 
as the Time=0 reference in Figure [2} and use this to calcu- 
late the burst time of arrival corrected to the Solar System 
Barycenter at infinite frequency, reported in Table For 
B2, since there are two clear burst components, we fit a 2- 
dimensional Gaussian to each and determined the Time=0 
reference as the center of the means of those Gaussians. 
Note that we do not fit pulse broadening functions to esti- 
mate the scattering timescale since the estimated scattering 
timescale from the Milky Way ISM is ~50ns (« 8 us) at 
1.4 GHz! and is consistent with the frequency structure 
we attribute to scintillation (see below). 


To more accurately determine the burst widths, we per- 
formed a 2-dimensional autocorrelation of the dynamic 
spectra, and fit these with 2-dimensional Gaussian functions 
(Extended Data Figure [2). Note the zero-lag noise spike is 
removed from this ACF. We convert the standard deviation 
of this Gaussian fit (in both time and frequency) to a full- 
width at half-maximum (FWHM) by multiplying the stan- 
dard deviation by the factor 2 V2In(2). We report the burst 
time width (tyiq = 1/ V2 x FWHM) and frequency extent 
(Ywia = FWHM) in Table In Figure [2| the colored bars 
below each burst profile indicate the +2o width used for 
calculations of the fluence and polarization fractions. 

In addition to the frequency extent measured with the 
2-dimensional ACF, there is another, narrower frequency 
scale evident in the one-dimensional frequency ACF (or 
simply by eye in the dynamic spectra of the bright bursts). 
Shown in Extended Data Figure [3A is the frequency ACF 
after subtracting the larger frequency scale Gaussian. We 
perform a least-squares fit of a Lorentzian function to the 
center of this ACF (defined by eye using the clearly visi- 
ble central feature in the ACF; Extended Data neuer 
A Lorentzian frequency ACF is expected for scintillation?! 
The fit function is of the form 


a 


+ b, (5) 


x? + Veit 
where a is the amplitude, b is a vertical offset, and Vecate 
is the scintillation bandwidth (defined as the half-width 
at half-maximum of the ACFE}, The scintillation band- 
width measurements are reported in Table [I] The Galactic 
ISM is expected to introduce a broadening due to scatter- 
ing of ~ 50ns (at 1.4 GHz) along this line-of-sight?! in 
rough agreement with our measured scintillation bandwidth 
(1/(27AyYgcint) ~ 27 ns). We therefore attribute this narrower 
spectral structure to scintillation from the Milky Way ISM. 
This interpretation is supported by the stronger correlation 
in the spectrum of two bursts (B1 and B2) separated by 4.3 
minutes, compared to the lack of correlation between bursts 
B3 and B4, which are separated by 2.5 hrs (Extended Data 
Figure [TO). The expected scintillation time is ~ 10 min- 
utes (at 1.4 GHz) at high Galactic latitudes, which is depen- 
dent on an effective velocity (the assumption for Galactic 
sources is 100 km s~! 28184] This effective velocity is likely 
smaller for an extragalactic source than for a Galactic pulsar 
in the same line of sight, since the effective velocity is usu- 
ally dominated by the velocity of the pulsar. Assuming an 
effective velocity of ~ 30 km s7}, the expected scintillation 
time will be ~ 30 minutes (at 1.4 GHz), consistent with the 
decorrelation timescale constraints using FRB 20200120E. 


Polarization calibration 


To study the polarization properties of the bursts, we must 
first calibrate the data. In these observations, we did 
not record a noise diode scan to use for delay calibra- 
tion; instead, we use the test pulsar data, also taken to 
inspect the general data quality, to calibrate the polari- 
metric data. This analysis strategy has been used suc- 
cessfully in previous work to determine the polarimet- 
ric properties of radio bursts from the Galactic magnetar 


SGR 1935+2154, detected using voltage data of the West- 
erbork RT1 VLBI backend#! and also in a study of the 
repeating FRB 20180916B, using data from the Effelsberg 
telescope as part of an EVN campaign! Since we detected 
bursts on separate epochs, we calibrate the data using the 
test pulsar scan closest in time to the bursts. B1 and B2 
are < 1hr from the PSR B0355+54 calibration scan, and 
B3 and B4 are ~ 4hr and < 2hr from the PSR B0355+54 
scan, respectively. We note that we could not recover the 
polarimetric properties of burst B5 likely due to the low S/N 
of the burst, and so the April 28 epoch is omitted from the 
remainder of this section. 

The data could exhibit leakage between the polariza- 
tion channels, which we assume only significantly affects 
Stokes V (defined as V=LL—RR using the PSR/IEEE con- 
vention), In the pulsar data, without calibrating the 
leakage, we reproduce the circular polarization fraction to 
within 3% and 1% of the published values for the Feb 
20 and Mar 7 epochs, respectively. We, therefore, apply no 
leakage calibration. 

A delay between the polarization hands is more crucial 
to correct for since it can significantly impact our inter- 
pretation of the linear polarization fraction and RM. Using 
the known RM of PSR B0355+54, 79 rad m2 we per- 
formed a brute-force search for the delay between polar- 
ization hands, D, that maximises the linear polarization, by 
rotating the data using the factor 

e72iRM(c* W) e72ivD, (6) 
where c is the speed of light and v is the frequency in Hz. 
We searched for delays between —15 and 15 ns, in steps of 
0.01 ns. Note that previous work with the VLBI recorder at 
Effelsberg showed an instrumental delay of 5.4 ns/6| The 
estimated delays are —0.18 ns and —4.11 ns for the Feb 20 
and Mar 7 observations, respectively. 

After removing the effect of the estimated delays, we per- 
formed RM synthesis on the burst data to estimate the 
RMs of the bursts. The Faraday spectra for each burst are 
shown in Extended Data Figure (itp. 

To refine the RM and delay measurements further, we 
performed a joint least squares fit of Stokes Q and U spectra 
normalised by the linear polarization L= -/Q? + U?, using 
the following equations 


QIL = cos(2(c7RM/v + vaD + 4)), (7) 


(8) 


where ¢ is a linear combination of the absolute angle of po- 
larisation on the sky (referenced to infinite frequency) and 
the phase difference between the R and L polarisation chan- 
nels. The estimates of the delay and RM are given as initial 
guesses to the fit, and we simultaneously fit the pulsar and 
the two bursts from the same observational epoch. We force 
the delay to be the same between the pulsar and FRB data, 
fix to the known RM of the pulsar (79 rad m°), and ini- 
tially we allow for a different RM per burst. The measured 
delays are —0.22 + 0.02 ns and —4.16 + 0.03 ns, for the Feb 
20 and Mar 7 observations, respectively, and the measured 
RMs are presented in Table The QU-fits are shown in 
Extended Data Figure The measured RMs have large 


U/L = sin(2(c7RM/v” + vaD + ¢)), 
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uncertainties due to covariances between the fit parameters, 
which is difficult to alleviate due to the low number of rota- 
tions (due to either RM or delay) across the burst frequency 
extent. We conclude that the bursts in this work have consis- 
tent RM values (B2’s RM is slightly over 30 from B3, but 
this apparent difference should be verified in future observa- 
tions of this source using an independent delay calibrator), 
which is also in agreement with the previously reported RM 
(—29.8 rad m24 We note that neither this work nor Bhard- 
waj et al.4! correct for the ionospheric contribution to the 
RM, but this effect is likely to be $ 2 rad m7. 

The data for each burst were then corrected for the mea- 
sured delay and RM. The polarization position angle (PPA) 
was calculated using the following equation: 


PPA = 0.5 arctan (5) i (9) 


and, to correct for the parallactic angle, this is rotated by 


gaspi sin(HA) cos(iat) 
7 (sin( iat) cos(ô) — cos(a) Sin(6) cos(HA) / 
(10) 


where HA is the hour angle of the burst, jat is the lati- 
tude of the Effelsberg telescope, and 6 is the declination 
of FRB 20200120E. Due to the fact that our observations 
did not feature an independent polarization calibrator scan, 
there remains an uncalibrated absolute phase offset in the 
data. Therefore, we cannot compare the PPAs between our 
two observational epochs, or with other PPA measurements 
of FRB 20200120E. By allowing for individual RM values 
per burst, we find that B1 and B2 exhibit a APPA of 28.5°, 
and the APPA between B3 and B4 is 21.7°. The unbiased 
linear polarization is computed following Everett & Weis- 
berg, where 


-1, ifs > 1.57 


otherwise 


(11) 


Lunbias = 


where Lmeas = yQ? + U2, and g; is the standard deviation 
in the off-burst Stokes I data. 

Bursts B1 — B4 are highly linearly polarized (> 90%), 
and exhibit little-to-no evidence for circular polarization. 
There is a tentative 3 — 4c detection of 13 % and 6 % circu- 
lar polarization in B2 and B4, respectively. The linear and 
circular fractions, as well as the PPA offset from a weighted 
mean PPA of bursts from a single epoch, are presented in 
Table[2| and the polarization profile and PPA for each burst 
is shown in Figure [2] 

It has been seen in the literature that some repeating 
FRBs exhibit a constant PPA (A(PPA) < 10°) and RM per 
observing epoch 4116] Therefore, in addition to the individ- 
ual burst RM measurements, we fit for a global RM per 
observing epoch. In this global fit, we find that the PPA of 
B1 and B2 differ by 68.3° (Table [2). This large APPA may 
result from the low S/N of burst B1. The APPA between B3 
and B4 in this global fit is 23.6°, comparable to the differ- 
ence we measured in the independent-RM fit above, unsur- 
prising given the very similar RM values in the independent 
fits. 
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Extended Data Figure 1: The probability that the brightest 1-2 bin features in the 31.25 ns burst profiles are consistent 
with the local amplitude modulated noise distribution. Panels a—c show the 31.25 ns resolution profiles (grey) with off 
burst noise (black) shown for comparison. The burst name and time resolution is shown in the top left corner and the 
frequency range averaged over to produce the burst profile in the top right corner of the panels. Panels d-f are zoomed-in 
profiles containing the highest S/N feature in the burst profile. The colored region represent the local region (time span 
+1.5625 us) used to determine the probability density function (pdf; panels g-i) and cumulative density function (cdf; 
panels j-l). Note that the feature at the center of the colored region is not added to the distribution, since this is the 
feature which we want to determine the significance of relative to the local distribution. An exponential distribution fit 
is overplotted (colored lines) on the pdf and cdf. The highest S/N feature is represented by the vertical dashed line on 
the cdf (where in the case of B3, there are two dashed lines since the feature is 2 bins wide), and the legend shows the 
probability (or 1-cdf) of these features. The horizontal dashed lines on panels d-f represent the 3 o levels for single-bin 
features, using this local distribution (also for 2-bin features in the case of B3). 


17 


Frequency Lag [MHz] 


-0.3 -0.2 -0.1 0.0 


Frequency Lag [MHz] 


-0.10 -0.05 0.00 0.05 0.10 -03 -0.2 -0.1 00 01 02 03 


Time Lag [ms] 


Frequency Lag [MHz] 


-0.2 -0.1 0.0 0.1 0.2 
Time Lag [ms] 


Extended Data Figure 2: Low time resolution 2D autocorrelation functions (ACFs) of the bursts detected from 
FRB 20200120E. Panels f-j show the 2D ACF with colored contours overplotted representing the 2D Gaussian fit 1,2,3 
and 4 o. The zero lag spike is not plotted. The ACF is computed using filterbank data generated with SFX] with time 
and frequency resolution of 8 us and 125 kHz, respectively. The data were dedispersed using a DM of 87.75 pecm™’. 
Panels a—e show the frequency-averaged time ACF, with the frequency-averaged Gaussian fit overplotted, and similarly 
panels k—o show the time-averaged frequency ACF, with the time-averaged Gaussian fit overplotted. The time and fre- 
quency scales characterised in this plot, arise from the burst temporal width and frequency extent. The colored lines 
coordinate with other figures in this work (e.g. Figure [2Jand Extended Data Figure[I Ip. 
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Extended Data Figure 3: Measurement of the scintillation bandwidth in the autocorrelation function (ACF) of each of 
the five bursts from FRB 20200120E. Sub-figure a shows the time-averaged frequency ACF from Extended Data Figure[2] 
after subtracting the Gaussian fit (black). The Lorentzian fit to the central component is shown by the colored line. The 
burst name is shown in the top left of each panel, and the measured scintillation bandwidth (defined as the half-width at 
half-maximum of the Lorentzian; Rickett!83) is shown in the top right of each panel. Sub-figure b shows a colormap of 
the expected scattering timescale at 1.4GHz as a function of Galactic longitude and latitude, from the NE2001 Galactic 
electron density model28l The sky positions of both FRB 20200120E and FRB 20180916B are shown by the black cross 
and white plus, respectively. 


Extended Data Table 1: High time resolution autocorrelation function and power spectra results. 


Burst Characteristic timescales [us] Red. y7* Power law index” A Bayesian Information Criterion? Goodness of fit p-value’ Outlier p-value® 
B2 28.5 40.2 2.0 1.85 + 0.04 —20.5 0.50 0.98 
B3 0.04 + 0.3, 1.1 40.4, 35.6 +0.5 2.4, 1.64, 1.3 1.46 + 0.05 -12.1 0.53 0.99 
B4 0.05 + 0.02, 64.1 + 1.5 2.7,1.9 2.04 + 0.05 -15.4 0.66 0.97 


* The reduced y” of the Lorenzian fit to the high time resolution ACF in Extended Data Figures 4f] 

P The power law index measured for the power spectra. 

€ Metric for model comparison (BIC for power law red noise model — BIC for power law plus Lorentzian). 
4 Goodness of fit of the best fit model. 

© p-value of the highest outlier in the residuals of the power spectrum. 
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Extended Data Figure 4: Dynamic spectrum of burst B3 from FRB 20200120E with time resolution 31.25 ns, temporal 
autocorrelation function (ACF) and power spectrum (PS). Panel b shows the dynamic spectrum in the form of tempo- 
ral profiles per subband. This data was generated with SFXC, and each subband has been coherently dedispersed to 
87.7527 pc cm™°. Panel a shows the frequency-averaged burst profile. Panel ¢ shows the average power spectrum (PS) of 
the four subbands containing significant burst structure in the top panel (grey), with a downsampled PS (factor 8) over- 
plotted in black. The purple and blue lines represent fits to the PS of a red noise power law plus white noise model and a 
power law/white noise plus Lorentzian model, respectively. Panels d and e below show the residuals (2 x D/M, for data D 
and model M) of both models, matching the colors above. The dashed lines represent the perfect case of D = M. Panel f 
shows the average temporal ACF of the same four subbands (top panel, orange). For comparison the off burst ACF is also 
shown (grey). The residual of the average ACF subtracted the noise ACF is shown in panel g below, with the green and 
purple Lorentzian fits to the ACF residuals highlighting two distinct temporal scales in the data. Panel h shows a zoom in 
on the ACF residuals highlighting a third temporal scale by the cyan Lorentzian fit. 
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Extended Data Figure 5: The same as Extended Data Figure f] for burst B2 from FRB 20200120E. Note the average 
ACF only shows one temporal scale (unlike the three seen for burst B3). Additionally, we only plot the red noise power 


law plus white noise model since any wide Lorentzian features are less apparent in this power spectrum, than the case of 
B3. 
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Extended Data Figure 6: The same as Extended Data Figures [4] and [f5] for burst B4 from FRB 20200120E. Note the 
average ACF shows two temporal scales (unlike the 3 seen for burst B3). Additionally, we only plot the red noise power 


law plus white noise model since any wide Lorentzian features are less apparent in this power spectrum, than the case of 
B3. 


22 


1302-1462 MHz 0.5 


F 
ao 5 0.4 
o 
5 g 0.3 
a 5 c 
2 0.2 
& 
£ o1 
0 8 
0.0 
0 100 200 300 400 500 


1254-1430 MHz 
0.6 


0.4 


0.2 


Correlation coefficient 


0.0 


1302-1478 MHz 0.4 


0.3 
0.2 


0.1 


Correlation coefficient 


0.0 


0 100 200 300 400 500 600 0 10 20 30 0 50 
Time [ys] Time lag [us] 


Extended Data Figure 7: Correlation coefficient between single time bin (1 us) spectra of bursts B2, B3 and B4 from 
FRB 20200120E, testing for consistency with the Scintillating Amplitude Modulated Polarized Shot Noise prediction 
(0.33 for 100 % polarized emission=3)., Panels a—c show the 1 us burst profile with the burst name and time resolution 
shown in the top left corner and the frequency range averaged over to produce the burst profile in the top right corner of 
each panel. Panels d—e show the correlation coefficient between single time bin spectra above a S/N threshold of 9 as a 
function of the time separation between bins. The color gradient indicates the geometric mean of the two time bins used 
to determine the correlation coefficient (darker color implying a higher geometric mean). Also plotted is a histogram of 
the correlation coefficients. The colored line represents the geometric mean S/N weighted correlation coefficient. 
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Extended Data Figure 8: High time resolution polarimetric profile and polarization position angle (PPA) for burst B3 
from FRB 20200120E. Panels a—c show the PPA as a function of time, with the orange line representing the weighted 
best-fit line to the PPA. Only the PPAs above a linear S/N threshold of 5 are plotted. Panels d-f show the polarimetric 
profile of the burst sampled at 125 ns, with Stokes I (black), unbiased linear polarization (Everett & Weisberg, red) and 
circular polarization (blue). The yellow and blue regions plotted on panels a and d represent the time ranges used for 
plotting panels b,e and panels c.f, respectively. This data was generated with SFXC, with 4 MHz channels and coherently 
(within subbands) and incoherently (between subbands) dedispersed to 87.7527 pccm™>. The frequency information was 
averaged for the frequency range 1254 — 1430 MHz (visually, the extent of the burst in frequency), which, in this data 
product, corresponds to averaging by a factor of 44. 


24 


© — DM 87.7527 40 


1450 


Peak S/N 


1400 


1350 


Frequency [MHz] 


1300 


87.740 87.745 87.750 87.755 87.760 87.765 —0.04 -0.02 0.00 0.02 0.04 0.06 
Dispersion Measure [pc cm73] Time [ms] 


Extended Data Figure 9: Constraining the dispersion measure (DM) using the short timescale structure in burst B3 
from FRB 20200120E. Using data products generated using SFXC4] with a time and frequency resolution of 500 ns and 
1 MHz, respectively, we coherently (within subbands) and incoherently (between subbands) dedisperse to a range of DMs, 
downsample in time by a factor of 2, and compute the peak S/N of the frequency-averaged profile (panel a). Also plotted 
in panel a is a Gaussian fit to the peak S/N as a function of DM, with the best-fit DM (DM = 87.7527+0.0003 pc cm~?) 
shown by the grey line. Panel b shows the burst profile and panel c the dynamic spectrum coherently and incoherently 
dedispersed to the best-fit DM. The time and frequency resolutions used for plotting are 1 us and 16 MHz, respectively. 
The red lines on the dynamic spectrum indicate the frequency extent averaged over to determine the peak S/N per DM, 
and to produce the burst profile. 
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Extended Data Figure 10: Comparing the time-averaged spectra between bursts detected close in time. Panel a shows 
the time-averaged spectra for B1 (black) and B2 (grey). Panel b shows the time-averaged spectra for B3 (grey) and B4 
(black). In the top right of each panel, we quote the time between the two bursts in the plot, ót, and the correlation 
coefficient of the burst spectra, R. 
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Extended Data Figure 11: Rotation measure (RM) determination for B1 — B4 from FRB 20200120E. Panel a shows 
the Faraday spectrum per burst (labelled), determined using RM synthesis®4) with the grey line showing the previous RM 
measurement of this source#l The grey-scale scatter points panels b—e show the Stokes Q (top) and Stokes U (bottom) 
spectra normalised by the linear polarization L= yQ? + U?, where the darker color represents a higher S/N. The colored 
lines show the best-fit QU-fitting result, where we fit for both the RM and the delay between the polarization channels 
(Methods). The color of the line indicates which burst is being plotted, using the legend in the top figure. 
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